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We introduce the Peierls substitution to a two-dimensional discrete-time quantum walk on a square 
lattice to examine the spreading dynamics and the coin-position entanglement in the presence of 
an artificial gauge field. We use the ratio of the magnetic flux through the unit cell to the flux 
quantum as a control parameter. For a given flux ratio, we obtain faster spreading for a small 
number of steps and the walker tends to be highly localized around the origin. Moreover, the 
spreading of the walk can be suppressed and decreased within a limited time interval for specific 
rational values of flux ratio. When the flux ratio is an irrational number, even for a large number 
of steps, the spreading exhibit diffusive behavior rather than the well-known ballistic one as in the 
classical random walk and there is a significant probability of finding the walker at the origin. We 
also analyze the coin-position entanglement and show that the asymptotic behavior vanishes when 
the flux ratio is different from zero and the coin-position entanglement become nearly maximal in a 
periodic manner in a long time range. 

PACS numbers: 03.67.Ac,03.67.Bg,03.67.Lx,03.65.Vf,42.60.-v 


I. INTRODUCTION 

A quantum computer was envisioned by Feynman as 
a device overcoming the difficulty of simulating quantum 
mechanical systems with classical computers [l|. Today, 
ultracold atomic systems are of great interest for imple¬ 
menting highly controllable analogues of quantum sys¬ 
tems under consideration ill. One way of creating 
and controlling such systems is trapping ultracold neu¬ 
tral atoms in periodic potentials of optical lattices formed 
by a laser. Since the atoms are neutral, external electric 
or magnetic fields, which are essential for quantum phe¬ 
nomena such as the quantum Hall effect and topological 
phases, have no effect on their trajectories. On the other 
hand, analogous effects can be implemented into optical 
lattices artificially to extend their simulation abilities. 
There are several proposals ill and experiments in 
focusing on the creation of tunable artihcial gauge fields 
for ultracold neutral atoms in optical lattices by using 
atom-light interaction. 

The discrete-time quantum walk (QW) was originally 
proposed by Aharonov et al. as a quantum counterpart 
of the classical random walk, where the QW leads to bal¬ 
listic spread of the walker rather than the diffusive one 
observed in the classical case M- Quantum walks are 
useful for developing new quantum algorithms m and 
they provide a model for universal quantum computa¬ 
tion (l2l| . They also supply a fertile framework for simu¬ 
lating other quantum systems such as topological phases 
[T^ . Anderson localization [l3, nonlinear 5-kicked 
quantum systems [l^ , the breakdown of an electric-field 
driven system [13. and the creation of entanglement in 
bipartite systems [T^. Quantum walks can be realized 
experimentally in various physical s yste ms including ul¬ 
tracold atoms in optical lattices [isl - l^ . 
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Recent experimental studies on optical lattices can al¬ 
low the realization of QWs under artificial gauge fields. 
For example, it has been experimentally shown that the 
effect of an electric field on a charged particle can be 
mimicked by QWs in a one-dimensional optical lattice 
[ 2 ^. Also, a recent proposal by Boada et al. utilizes 
photonic circuits for the realization of QWs under an ar¬ 
tificial magnetic field [ 2 ^. In this paper we investigate 
the dynamics of two-dimensional (2D) QWs on a square 
lattice in the presence of an artificial magnetic field. For 
this purpose, we introduce position- and direction- de¬ 
pendent phases corresponding to the Peierls phases of the 
hopping terms between neighboring sites in the Hamil¬ 
tonian representing a charged particle under a uniform 
magnetic field. We can control the propagation of the 
walker by changing the magnetic flux $ through the unit 
cell. Depending on $, we show that ballistic behavior can 
be suppressed within a time interval or can be completely 
broken. It is known that, classical diffusive behavior in 
QWs is observed when quantum coherence is removed in 
some way, e.g., by decohering the coin and/or the posi¬ 
tion [23 ■ We show that QW is also diffusive at long times 
if the magnetic flux ratio a = <I>/$o (‘f’o being the flux 
quantum) is an irrational number and the walker remains 
highly localized at the origin throughout the walk. More¬ 
over, we demonstrate that when a is chosen properly the 
walk stops to propagate and propagates back towards 
the origin during a limited time interval. We also ana¬ 
lyze the entanglement between the coin and the position 
of the walker and show that the well-known asymptotic 
behavior vanishes when a 0. We observe that the coin 
and the position become maximally entangled at specific 
steps under the effect of the artificial magnetic field on a 
long time scale. 

This article is organized as follows. In Sec. |ll]we give 
a brief overview of QWs and introduce the Peierls model 
to the formalism. In Sec. imi we compare the behav¬ 
ior of the QW under rational and irrational a’s and we 
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demonstrate the localization of the walker. In Sec. llVI we 
examine the effect of the magnetic field on coin-position 
entanglement. In Sec. 0we summarize our results. 

II. PEIERLS MODEL IN QWS 

In analogy to the classical random walk, the master 
equation for the discrete-time quantum walk is given by 
where W is a unitary transformation 
describing the time evolution of the state vector j'l't) in 
discrete bipartite coin-position Hilbert space ® 
spanned by {|0)^, |1)^} and {\n)^ \ n S Z}, respectively. 
The operator W = S{C®X) is called a step of the walk and 
it is composed of a shift operator S and a coin operator 
C. The coin operator acts only on the coin space and 
it can be any unitary operation in SU(2). In the first 
proposal of the QW [l^ , C was chosen as the Hadamard 
gate 

Ch = ^(I0)c(0| + |0),(I| + |1),(0| - |I),(I|), (1) 

which is the one we use throughout the paper. The con¬ 
ditional shift operator on a line is 

+ 00 

(|0)c(0| 0 |n-l)p(n|-k |l)^(l| (8) |n-M)p(n|), 

n— — oo 

( 2 ) 

which acts on both spaces (coin and position) and it 
moves the walker to the left (right) when the coin compo¬ 
nent is in the state |0)^ (|l)c)' If walk starts with the 
initial state I'I'o); after t steps, the final state becomes 

|Tt) =W‘|^'o) = ^a„.e;t|c,n) (3) 

c,n 

where an,c;t are the site amplitudes and c S {0,1}. 
The probability of being found at any position Pn-,t = 
|aTi,c;tP is calculated by summing over the probabili¬ 
ties in the coin space. 

Quantum walks on a line can be extended to higher di¬ 
mensions by enlarging the coin and the position spaces. 
For 2D QWs, the coin space can be chosen as four di¬ 
mensional with the basis {|0)^, |2)^, |3)^}, while the 

position space is spanned by {|n,m)p \ n,m € Z}. The 
coin can be interpreted as a single four-level coin or two 
different two-level coins. For both cases, the coin op¬ 
erator can be factorizable, i.e., C = Ch <8) Ch, or non- 
factorizable, e.g., the Grover coin (the so-called Grover 
walk) or the discrete Fourier transform coin [2^. When 
the walker component |n, to)^ corresponds to the position 
eigenstates of the walker on a square lattice, the shift op¬ 
erator is defined as a single operation that moves the 
walker in the left down, left up, right down and right up 
directions (i.e., towards corners) for respective coin states 



FIG. 1. Phases gained by the walker while hopping between 
neighboring sites in the y direction. These phases are anal¬ 
ogous to those acquired by a charged particle subjected to 
a constant external magnetic field in the 2 direction, in the 
Landau gauge. Black circles denote the sites of the square lat¬ 
tice labeled by the integer pairs (n, m). The arrows show two 
possible trajectories. The walker gains a total phase of 
along the smaller paths and along the larger paths. 

In general, for an arbitrary closed path, the coefficient of a 
gives the number of unit cells enclosed. No phase is gained 
along the x direction. 

[^ . On the other hand, experimentally, the walker can 
step to the nearest neighbors rather than the corners. 
Therefore, if we want the walker to be found at the cor¬ 
ners in one step, we have to implement the shift operator 
as two separate operations, i.e., a shift along one axis 
followed by a second shift along the other axis. An alter¬ 
native method for realization of 2D QWs is to use a single 
two-level coin instead of a four-level coin (sain. In this 
scheme, a step is defined asU = Sy{C®X)Sx{C®I), where 
the walker is first shifted along the x direction followed by 
a shift along y direction. The Hadamard operator in Eq. 
(ED is chosen as the coin operator and it is applied before 
each shift. It has been shown that, with this scheme, 
the probability distribution of the Grover walk can be 
mimicked. This alternate scheme has advantages over a 
walk with a four-level coin when experimental aspects for 
square lattices are considered. Therefore, in this paper, 
we use this alternate scheme. 

We consider QWs on a square optical lattice with a 
tunable artificial gauge field. We label the space coor¬ 
dinates as X = na and y = ma, where a is the lat¬ 
tice constant. We choose the symmetrical initial state 
l^o) = l^)c ® where both n and m are defined as 

zero and |'(/')c = "^( 10)0 + *|I)c)- introduce the shift 
operator in the y direction as 

Sy =|0)^(0| i8)^e+*^™”|n,m- I)p(n,m| 

( 4 ) 

n,7n 
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(a) (b) 



FIG. 2. (a) Spreading of the walk after a different number of steps t with respect to the magnetic flux <1? in units of flux 
quantum $o. The symmetrical initial state |^o) is used. Only one period of variances is drawn within the interval [0, 1]. The 
vertical axis is rescaled with the maximum value of the variance for each t. For the rational values of a, sharp peaks become 
apparent when the number of steps is increased, (b) Spreading of the walk with respect to the number of steps for different 
values of a. When a is the golden ratio f#, we used two different initial conditions. The lower solid line corresponds to the 
symmetric case used for the other a values and the upper one is |^o) = |0}^ ® |0, 0)^. 


where a £ [0,1] is the tuning parameter and 
are both site and direction-dependent phases of hopping 
terms between neighboring lattice sites. At this point 
our approach differs from [25^ and [^ . where only site- 
dependent phases on a one-dimensional lattice are used. 
When we consider the motion of a charged particle on a 
square lattice under a uniform magnetic field B = B^z, 
the appropriate Peierls substitution is given by 

t—J A.dl^ (5) 

where t is the nearest-neighbor hopping amplitude and 
$0 = ~ 4.14 X 10“^® Wb is the flux quantum {h 

and e being the Planck constant and elementary charge, 
respectively). Here ri and t -2 denote the initial and final 
positions of the particle, respectively, and the integral is 
evaluated along the line connecting these points. If the 
vector potential is chosen as A = B^xy (Landau gauge), 
the transition amplitude along the x direction remains 
unaffected while along the y direction it gains a factor of 

exp ^±i27r—(6) 

when hopping from site m to mil. Here $ = Bqo?' 
is the magnetic flux through the unit cell. Therefore, 
we interpret the phases in Eq. o as artificial Peierls 
phases shown in Fig. [TJ In this paper, we consider both 
rational and irrational flux ratios a = 4>/<i>o to find out 
the effects of the magnetic field on the characteristics of 
QWs such as the variance, participation ratio, and coin¬ 
position entanglement. 


III. RATIONAL VS IRRATIONAL FLUX 
RATIOS 

We consider the variance tr^ of the probability distri¬ 
bution, 

+ 'rn^)P{n,m)-t. (7) 

n^m 

as a measure of the spreading. Here, P(n,m)-,t is the prob¬ 
ability of the walker being found at site (n, m) after t 
steps and we look for how the variance behaves under 
rational and irrational magnetic fields. When the flux 
ratio a is an integer, the original translational symme¬ 
try of the lattice is preserved. Similarly, when it is a 
rational number such as p/?, where both p and q are co¬ 
prime integers, translational symmetry is preserved only 
if the unit cell is considered q times as large. However, 
for irrational a values, the number of unit cells enclosed 
by the walker (see Fig. [1} is incommensurable with the 
parameter a. Therefore, we cannot exploit a rescaling 
as we did in the rational case. In Fig. [H[a), we show 
the change of the variance with respect to a for different 
step numbers. Although the variance is meaningful for 
the walks with a large number of steps, we intentionally 
present the results for a small number of steps to demon¬ 
strate an interesting effect of magnetic fields on QWs. 
Note that since the variance repeats itself with a period 
of 1/2 within the interval [0,1], only one period is drawn 
for each case. When t = 2,4, or 8, for each case there are 
two maxima at non-zero a values. While the number of 
steps is increasing, the maxima move to the left and their 
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(a) 


(b) 




Number of Steps (t) 


FIG. 3. (a) Probability distribution of the walk after 100 steps for different flnx ratios. Here is the golden ratio, (b) 
Probability of finding the walker in the vicinity of the origin after t steps: the snm of the probabilities at sites where 

n,m (z {—2, 0, 2} with respect to the nnmber of steps. Since an odd number of steps gives vanishing probability, the probabilities 
for only an even nnmber of steps are shown. Data are smoothed out via Gaussian filtering for simplicity. The inset shows a 
comparison of the actual data with the filtered data for a = ^. 


positions converge to a = 0 and a = 1/2. For example, 
the analytic expression of the variance for two steps is 

0-2 = 3 + 2cos^(27ra — ^). (8) 

It is clear that the maximum value of Eq. (jS]) occurs at 
a = 1/8 and 5/8. It is also notable that in this case the 
probability at the origin, 

-P(0.0);2 = ^ cos^(27rQ; + ^), (9) 

becomes zero because of the destructive interference of 
the incoming amplitudes. In other words, the walker 
avoids stepping into the center at the second step as a 
result of the applied field. Similarly, maxima for four 
and eight steps occur at a = 3/100 and a = 3/400, re¬ 
spectively. Hence, with an appropriate choice of the a, 
the walk spreads faster than the case where there is no 
applied field for a small number of steps. For 20 and 60 
steps, spreading of the walk is dramatically reduced by 
the non zero values of a except 1/2. However, some peaks 
occur at the rational values of a and more peaks become 
apparent when the number of steps is larger. This hap¬ 
pens because the more steps taken, the more routes that 
satisfy translational symmetry corresponding to different 
a values can be followed. Translational symmetry due to 
any rational number a requires the walker to take at 
least a number of steps that is enough to follow a closed 
path covering the unit cell rescaled by a appropriately. 
Therefore, we increase the a resolution of the walk by 
increasing the step number and thus the walk becomes 
more sensitive to the a values. By counting the number 
of peaks in Fig. [^Ja) roughly, we see that while the sys¬ 
tem is able to resolve 4 of those symmetries at 20 steps, it 


can resolve more than 20 if we increase the step number 
to 60. 

In Fig. H a) we see that if a = 1/2, the variance is 
the same as the case where there is no magnetic flux 
a = 0 for a large number of steps. Actually, the walk 
has exactly the same dynamics for this two case and the 
same probability distribution is obtained after each step. 
The reason is that in the second and each succeeding 
step, translational symmetry is preserved, i.e., in each 
step the gained phases cancel each other after following 
a closed path that covers two unit cells. Hence, there 
is no net effect of magnetic flux ol = 1/2 on the walk. 
Note that no components of the wave function interfere 
with each other at the first step, i.e., the wave function 
only spreads towards the nearest corners. Therefore, the 
possible effects of the magnetic field occur after the first 
step for each walk under our consideration. 

In Fig. [2jb) we demonstrate that, for two different 
initial states, the walk spreads diffusively in a long time 
range if a is an irrational number and the golden ratio 
= ^~^ ), due to the broken translational symmetry 
of the lattice. In this case, there is no rescaled unit cell 
that is commensurable with the original one to obtain a 
translational symmetry. In contrast, if a is rational, after 
few steps, spreading of the walk returns to its usual ballis¬ 
tic behavior. A remarkable case in which the spreading 
is nearly stopped and reversed during steps 35-50 and 
50-85, respectively, occurs for a = 21/44. Also, further 
simulations show that similar effects as in a = 21/44 can 
be observed if a is chosen sufficiently close to 0 or 1/2. 
In general, we can conjecture that when a is a rational 
number, even if the spreading of the walk fluctuates at 
the beginning, it will revert to its original ballistic be¬ 
havior after a finite number of steps. Only irrational a 
values permanently suppress the spreading and result in 
a diffusive behavior. 
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(b) 



Number of steps (t) 


FIG. 4. (a) Asymptotic values of coin-position entanglement of a QW on a square lattice with two-level initial coin states 
= cos ||0) -I- e“^sin ||1) for a — 0. The figure is symmetrical about rj) = tt and varies between its maximum ~ 0.903 and 
minimum ~ 0.945. The maximum value is attained for symmetrical initial coin states = ( 7 r/ 2 , 7 r/ 2 ) and ( 7 r/ 2 , 37 r/ 2 ). (b) 

Time dependence of the entanglement for specific values of the initial coin state. The asymptotic behavior can be observed 
explicitly. The insets show that when q 7 ^ 0, the asymptotic behavior vanishes. For a = (left inset) and for q = 1/4 (right 
inset), after the fourth step, the coin and the position become almost maximally entangled {Se — 0.99) every two steps. 


Figure [3Ka) demonstrates that whether a is a rational 
or an irrational number, the walker tends to be local¬ 
ized around the origin. However, some values of a cause 
stronger localization. In Fig. [3Kb) the sum of probabil¬ 
ities of the walker being found at the origin and some 
nearby sites with respect to the number of steps is given. 
When a is a rational number, even though the robustness 
of the localization against step number changes for dif¬ 
ferent a’s, the probability of finding the particle around 
the origin converges to zero while the number of steps 
is increasing. In contrast, when a = the probability 
does not converge to zero within the interval [0,1000] 
and moreover it increases to approximately 0.3 when 
t = 1000. Although we examine the sum of probabili¬ 
ties at several sites around the origin. Fig. |3Ka) ensures 
that the significant contribution to this sum comes from 
the origin. Therefore, an irrational flux ratio guarantees 
that the probability is the highest at the origin in the 
range up to 1000 steps. 

Although our results include only one irrational num¬ 
ber (the golden ratio), further simulations for I/tt, 1 /e, 
l/-\/2, and 1/C(3) show that all of them exhibit diffu¬ 
sive behavior on average. As we have mentioned above, 
the reason is the incommensurability of the number of 
unit cells enclosed by the walker and the magnetic flux 
ratio for each case. Note that, as shown in Fig. HKb), 
the walk can exhibit higher spreading rates temporar¬ 
ily, e.g., between t = 250 and 600 for a = ^, but on 
average, spreading fluctuates around the diffusive trend. 
However, these temporary deviations from the diffusive 
spreading can extend over relatively large time intervals 
for some irrational numbers. Among our simulations, the 
only example is a = I/tt, where we observe such a devi¬ 


ation over the interval from t = 100 to 2000. A possible 
explanation of this behavior can be given by answering 
the question of how well an irrational number can be ap¬ 
proximated by the rationals. It is known that the best ra¬ 
tional approximations to an irrational number are found 
in its convergents of continued fraction which are repre¬ 
sented by Ci = [oo.ai,--- ,ai] = pt/q^ [1^. By adding 
more terms, we obtain better approximations. Since the 
Peierls phases can be written as where 

tterr,! = |ci — a|, the phases allow various transla¬ 

tional symmetries determined by qi. We also know that 
if aerr,i IS Sufficiently close to 0 or 1/2 , the effects re¬ 
sult from cannot be observed for a while (see 

a = 21/44 in Fig. |3|and a = 1/500 in Fig. [3]). While the 
number of steps increases, the walk continuously tries to 
catch different translational symmetries required by the 
convergents but due to the accumulation of 
phases, this can never be done completely, which thereby 
results in temporary higher spreading rates. Therefore, 
if the Oerr’s are sufficiently small for a given irrational a, 
we can expect to observe higher spreading rates extend¬ 
ing over relatively long time intervals. What makes I/tt 
special is that its convergents have the minimal errors 
when we compare it with the other irrational numbers 
we used. It has the best rational approximations among 
the others and this can be the reason for the deviations in 
long ranges. A detailed analysis including a careful com¬ 
parison of the convergents of different irrational numbers 
and the corresponding errors is left to future research. 
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IV. EFFECT OF MAGNETIC FIELDS ON 
ENTANGLEMENT 

The shift operator generates entanglement between the 
coin and the position degrees of freedom [Sdl - fs^ . Since 
the step operator U is unitary, the density matrix pt = 
at any step t will be a pure state. Therefore, 
it is convenient to use the von Neumann entropy as a 
measure of the coin-position entanglement 

=-^AJog2 Ai, (10) 

i 

where the Ai’s are the eigenvalues of the reduced density 
matrix pc = tr^ (pt) obtained by tracing out the position 
degree of freedom. The asymptotic behavior of the entan¬ 
glement in QWs on a line and 2D QWs with a four-level 
coin is already known [s^, . Here we demonstrate that 

a 2D QW with a two-level coin also exhibits the same be¬ 
havior for all initial coin states for a = 0. In Fig. H^a) the 
dependence of the asymptotic values of coin-position en¬ 
tanglement on initial coin states is given. It can be seen 
that entanglement varies between a maximum of 0.945 
and a minimum of 0.903. The maximum value occurs for 
the symmetrical initial state |dtg). In Fig. HKb) time de¬ 
pendence of coin-position entanglement for three specific 
initial coin states is given. For all cases, entanglement 
decays to an asymptotic value. 

The inset in Fig. Hb) shows the change in coin¬ 
position entanglement in the presence of a magnetic field 
(for a ^ 0 and a ^ 1/2). For both a = 1/4 and 
a = ‘S^ the asymptotic behavior disappears completely. 
Although we show only two cases for simplicity, we ob¬ 
serve similar behaviors for the other a values. When 
a= 1/4, the value Se — 0.99 indicates that the coin and 
the position are almost maximally entangled in every odd 
step after the fourth step. When a = i#, we obtain al¬ 
most maximally entangled states in a quasi-periodic man¬ 
ner in a long time range again. To gain insight into the 
cause of the large values of entanglement, we consider the 
participation ratio at step number t, 

\n,m / 

which can be interpreted as an estimate for the number of 
sites over which the walker is distributed. While V = 1 
indicates that the walker is completely localized at only 
one site, N = d indicates a uniform spreading over d sites. 
In Fig. [5]we compare the coin-position entanglement and 
the participation ratio for a small magnetic field. During 
the first 10 steps, we see that the entanglement is exactly 
the same for a = 0 and a = 1/500. Further simulations 
for smaller values of a show that the entanglement is the 
same as the a = 0 case even for a larger number of steps. 
This observation suggests that the effects of the magnetic 
field start to appear after a sufficient amount of phases 



FIG. 5. Comparison between coin-position entanglement 
(thick) Se and participation ratio (thin) N when a = 1/500. 
Dashed lines represent the behavior of each case when q = 0. 

has accumulated. After about 400 steps, the participa¬ 
tion ratio starts to oscillate strongly and we also observe 
quasi periodic oscillations in Se-, which results in almost 
maximally entangled states as in the insets of Fig. HKb). 
Since strong oscillations in the participation ratio imply 
the tendency of the walker to return to the origin, we can 
conclude that localization of the walker plays an impor¬ 
tant role in the large values of coin-position entanglement 
in the presence of a magnetic field. 


V. CONCLUSION 

We have studied the spreading properties and coin¬ 
position entanglement for the QW on a square lattice 
under an artificial magnetic field. We have demonstrated 
that the presence of such fields increases the spreading 
for a small number of steps. Moreover, irrational flux 
ratios cause a diffusive spreading rather than a ballistic 
one even in the long time range because of the broken 
translational symmetry. This also causes the probability 
at the origin to be the highest even for a large number 
of steps. For rational flux ratios, ballistic spreading can 
be suppressed within a limited time range. However, the 
walk returns to the original ballistic behavior after a finite 
number of steps. 

We have demonstrated that the coin-position entangle¬ 
ment in a QW on a square lattice with a single two-level 
coin exhibits an asymptotic behavior as in the four-level 
case. We have also shown that this behavior changes in 
the presence of an artificial gauge field and it is possible 
to keep the coin and the position maximally entangled 
in a long time range if the field is chosen appropriately. 
However, we have not observed any distinguishing prop¬ 
erty of irrational or rational magnetic flux ratios while 
examining the coin-position entanglement. 

Finally, we note that our scheme may be realized with 
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today’s technology and some of our results, especially 
the ones that require only a few number of steps, may be 
verified. We believe that our work can provide a further 
step towards simulating many-body quantum systems in 
gauge fields by engineering the interactions of ultracold 


atoms with light. 
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